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Lecture 1 



Nonlinear optics is the discipline in physics in which the electric polarization density of the medium 
is studied as a nonlinear function of the electromagnetic field of the light. Being a wide field of 
research in electromagnetic wave propagation, nonlinear interaction between light and matter leads 
to a wide spectrum of phenomena, such as optical frequency conversion, optical solitons, phase 
conjugation, and Raman scattering. In addition, many of the analytical tools applied in nonlinear 
optics are of general character, such as the perturbative techniques and symmetry considerations, 
and can equally well be applied in other disciplines in nonlinear dynamics. 

The contents of this course 

This course is intended as an introduction to the wide field of phenomena encountered in nonlinear 
optics. The course covers: 

• The theoretical foundation of nonlinear interaction between light and matter. 

• Perturbation analysis of nonlinear interaction between light and matter. 

• The Bloch equation and its interpretation. 

• Basics of soliton theory and the inverse scattering transform. 

It should be emphasized that the course docs not cover state-of-the-art material constants 
of nonlinear optical materials, etc. but rather focus on the theoretical foundations and ideas of 
nonlinear optical interactions between light and matter. 

A central analytical technique in this course is the perurbation analysis, with its foundation 
in the analytical mechanics. This technique will in the course mainly be applied to the quantum- 
mechanical description of interaction between light and matter, but is central in a wide field of 
cross-disciplinary physics as well. In order to give an introduction to the analytical theory of 
nonlinear systems, we will therefore start with the analysis of the nonlinear equations of motion 
for the mechanical pendulum. 

Examples of applications of nonlinear optics 

Some important applications in nonlinear optics: 

• Optical parametric amplification (OPA) and oscillation (OPO), Tiuj-p — > TwOg + Ty^i- 

• Second harmonic generation (SHG), Hlo + Hlo ^ h{2Lo). 

• Third harmonic generation (THG), hw + fiw + fiu) ^ Ji{3u}). 

• Pockels effect, or the linear electro-optical effect (applications for optical switching). 

• Optical bistability (optical logics). 

• Optical solitons (ultra long-haul communication). 
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A brief history of nonlinear optics 

Some important advances in nonlinear optics: 

• Townes et al. (1960), invention of the laser. ^ 

• Franken et al. (1961), First observation ever of nonlinear optical effects, second harmonic 
generation (SHG).^ 

• Terhune et al. (1962), First observation of third harmonic generation (THG).^ 

• E. J. Woodbury and W. K. Ng (1962), first demonstration of stimulated Raman scattering.* 

• Armstrong et al. (1962), formulation of general permutation symmetry relations in nonlinear 
optics.^ 

• A. Hasegawa and F. Tappert (1973), first theoretical prediction of soliton generation in optical 
fibers.^ 

• H. M. Gibbs et al. (1976), first demonstration and explaination of optical bist ability.'' 

• L. F. Mollenauer et al. (1980), first confirmation of soliton generation in optical fibers.^ 

Recently, many advances in nonlinear optics has been made, with a lot of efforts with fields of, 
for example, Bose-Einstcin condensation and laser cooling; these fields arc, however, a bit out of 
focus from the subjects of this course, which can be said to be an introduction to the 1960s and 
1970s advances in nonlinear optics. It should also be emphasized that many of the effects observed 
in nonlinear optics, such as the Raman scattering, were observed much earlier in the microwave 
range. 

Outline for calculations of polarization densities 

Metals and plasmas 

From an all-classical point-of-view, the calculation of the electric polarization density of metals and 
plasmas, containing a free electron gas, can be performed using the model of free charges acting 
under the Lorenz force of an electromagnetic field, 

Cp'Y dv 

where E and B are all-classical electric and magnetic fields of the electromagnetic field of the light. 
In forming the equation for the motion of the electron, the origin was chosen to coincide with the 

center of the nucleus. 

^ Charles H. Townes was in 1964 awarded with the Nobel Prize for the invention of the ammonia 

laser. 

^ Franken et al. detected ulvtraviolet light (A = 347.1 nm) at twice the frequency of a ruby laser 
beam (A = 694.2 nm) when this beam traversed a quartz crystal; P. A. Franken, A. E. Hill, C. W. 
Peters, G. Weinreich, Phys. Rev. Lett. 7. 118 (1961). Second harmonic generation is also the first 
nonlinear effect ever observed where a coherent input generates a coherent output. 

^ In their experiment, Terhune et al. detected only about a thousand THG photons per pulse, 
at A = 231.3 nm, corresponding to a conversion of one photon out of about 10^^ photons at the 
fundamental wavelength at A = 693.9 nm; R. W. Terhune, P. D. Maker, and C. M. Savage, Phys. 
Rev. Lett. 8, 404 (1962). 

^ E. J. Woodbury and W. K. Ng, Proc. IRE 50, 2347 (1962). 

^ The general permutation symmetry relations of higher-order susceptibilities were published by 
J. A. Armstrong, N. Bloembergen, J. Ducuing, and P. S. Pershan, Phys. Rev. 127, 1918 (1962). 

^ A. Hasegawa and F. Tappert, "Transmission of stationary nonliner optica pulses in dispersive 
optical fibers: I, Anomalous dispersion; II Normal dispersion", Appl. Phys. Lett. 23, 142-144 
and 171-172 (August 1 and 15, 1973). 

H. M. Gibbs, S. M. McCaU, and T. N. C. Venkatesan, Phys. Rev. Lett. 36, 1135 (1976). 

^ L. F. Mollenauer, R. H. Stolen, and J. P. Gordon, "Experimental observation of picosecond 
pulse narrowing and solitons in optical fibers", Phys. Rev. Lett. 45, 1095-1098 (September 29, 
1980) ; the first reported observation of solitons was though made in 1834 by John Scott Russell, a 
Scottish scientist and later famous Victorian engineer and shipbuilder, while studying water waves 
in the Glasgow-Edinburgh channel. 
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Dielectrics 

A very useful model used by Drude and Lorentz^ to calculate the linear electric polarization of the 
medium describes the electrons as harmonically bound particles. 

For dielectrics in the nonlinear optical regime, as being the focus of our attention in this course, 
the calculation of the electric polarization density is instead performed using a nonlinear spring 
model of the bound charges, here quoted for one-dimensional motion as 

+ + a^'^x, + a^^^xl + a^^^xl + ... = -eE,{t). 

As in the previous case of metals and plasmas, in forming the equation for the motion of the 
electron, the origin was also here chosen to coincide with the center of the nucleus. 

This classical mechanical model will later in this lecture be applied to the derivation of the 

second-order nonlinear polarzation density of the medium. 

Introduction to nonlinear dynamical systems 

In this section we will, as a preamble to later analysis of quantum-mechanical systems, apply 
perturbation analysis to a simple mechanical system. Among the simplest nonlinear dynamical 
systems is the pendulum, for which the total mechanical energy of the system, considering the 
point of suspension as defining the level of zero potential energy, is given as the sum of the kinetic 
and potential energy as 

E = T + V=^\p\^-mgl{cos^-l), (1) 

where m is the mass, g the gravitation constant, I the length, and •& the angle of deflection of the 
pendulum, and where p is the momentum of the point mass. 




Figure 1. The mechanical pendulum. 

From the total mechanical energy (1) of the system, the equations of motion for the point mass is 
hence given by Lagranges equations, 

d f dL\ dL ^ . 

where qj are the generalized coordinates, pj = (jj are the components of the generalized momentum, 
and L = T — V is the Lagrangian of the mechanical system. In spherical coordinates (p, ip, "&), the 
momentum for the mass is given as its mass times the velocity, 

p = ml{'&e^ + ipsin'detp) 

and the Lagrangian for the pendulum is hence given as 

^ = ^ (4 + P^) + mgl{cos 1? - 1) 
■ 

= ^(^^ + tp^sin^i?) + mgl{cos'd - 1). 



^ R. Becker, Elektronen Theorie, (Teubner, Leipzig, 1933). 

^° Herbert Goldstein, Classical Mechanics, 2nd ed. (Addison- Wesley, Massachusetts, 1980). 
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As the Lagrangian for the pendulum is inserted into Eq. (2), the resulting equations of motion are 
for qj = if and 'd obtained as 



+ 



~di 



sin ip cos = 0, 



and 



df^ 



+ {g/l)sm'd = 0, 



(3a) 



(3b) 



respectively, where we may notice that the motion in e,p and directions are decoupled. We may 

also notice that the equation of motion for ip is independent of any of the physical parameters 
involved in te Lagrangian, and the evolution of ip{t) in time is entirely determined by the initial 
conditions at some time t = to- 

In the following discussion, the focus will be on the properties of the motion of i^it). The 
equation of motion for the coordinate is here described by the so-called Sine-Gordon equation. 
This nonlinear differential equation is hardf to solve analytically, but if the nonlinear term is 
expanded as a Taylor series around = 0, 



dt^ 



l?3 



^5 

5[ 



0, 



Before proceeding further with the properties of the solutions to the approximative Sine-Gordon, 
including various orders of nonlinearities, the general properties will now be illustrated. In order to 
illustrate the behaviour of the Sine-Gordon equation, we may normalize it by using the normalized 
time r = [g/iy^'^t, giving the Sine-Gordon equation in the normalized form 



dr'^ 



sin I? = 0. 



(4) 



The numerical solutions to the normalized Sine-Gordon equation are in Fig. 2 shown for initial 
conditions (a) 2/(0) = 0.1, (b) y(0) = 2.1, and (c) y(0) = 3.1, all cases with j/'(0) = 0. 



m. v(Q)=0.1ll, Oy((lVdt=CJ 




m, y(0)=3.10, d¥(0)Wl=D 




Figure 2. Numerical solutions to the normalized Sine-Gordon equation. 



The term "Sine-Gordon equation" has its origin as an allegory over the similarity between the 
time-dependent (Sine-Gordon) equation 

d'^ip 1 d'^ip o . 

appearing in, for example, relativistic field theories, as compared to the time-dependent Klein- 
Gordon equation, which takes the form 

d'^ip 1 d'^ifi _ 2 



The Sine-Gordon equation is sometimes also called "pendulum equation" in the terminology of 
classical mechanics, 
t Impossible? 
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First of all, wc may consider the linear case, for which the approximation sin?? w holds. For 
this case, the Sine-Gordon equation (4) hence reduces to the one-dimensional linear wave-equation, 
with solutions i? = Asm{{g /lY/'^{t — to)). As seen in the frequency domain, this solution gives a 
delta peak at cj = [g/l)^^'^ in the power spectrum |?9(cl')P, with no other frequency components 
present. However, if we include the nonlinearities, the previous sine-wave solution will tend to 
flatten at the peaks, as well as increase in period, and this changes the power spectrum to be 
broadened as well as flattened out. In other words, the solution to the Sine-Gordon give rise to 
a wide spectrum of frequencies, as compared to the delta peaks of the solutions to the linearized, 
approximative Sine-Gordon equation. 

From the numerical solutions, we may draw the conclusion that whenever higher order nonlinear 
restoring forces come into play, even such a simple mechanical system as the pendulum will carry 
frequency components at a set of frequencies differing from the single frequency given by the 
linearized model of motion. 

More generally, hiding the fact that for this particular case the restoring force is a simple sine 
function, the equation of motion for the pendulum can be written as 

^ + + a(i)i? + a(2),?2 ^ ^(3)^3 + . . . = q. (5) 

This equation of motion may be compared with the nonlinear wave equation for the electromagnetic 
field of a travelling optical wave of angular frequency w, of the form 



which clearly shows the similarity between the nonlinear wave propagation and the motion of the 
nonlinear pendulum. 

Having solved the particular problem of the nonlinear pendulum, we may ask ourselves if the 

equations of motion may be altered in some way in order to give insight in other areas of nonlinear 
physics as well. For example, the series (5) that define the feedback that tend to restore the 
mechanical pendulum to its rest position clearly defines equations of motion that conserve the 
total energy of the mechanical system. This, however, in generally not true for an arbitrary series 
of terms of various power for the restoring force. As we will later on see, in nonlinear optics 
we generally have a complex, though in many cases most predictable, transfer of energy between 
modes of different frequencies and directions of propagation. 

The anharmonic oscillator 

Among the simplest models of interaction between light and matter is the all-classical one-electron 
oscillator, consisting of a negatively charged particle (electron) with mass rrie, mutually interacting 
with a positively charged particle (proton) with mass rUp, through attractive Coulomb forces. 

E0t5:=i^e;:::^:: " r:::::^:::::: 

' Z Z Z Z Z Z * Z 3~~c ^ — ^(*^c -^n) ' z z z z 

^^^^ 

— 




Figure 5. Setup of the one-dimensional undamped spring model. 

In the one-electron oscillator model, several levels of approximations may be applied to the 
problem, with increasing algebraic complexity. At the first level of approximation, the proton 
is assumed to be fixed in space, with the electron free to oscillate around the proton. Quite 
generally, at least within the scope of linear optics, the restoring spring force which confines the 
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electron can be assumed to be linear with the displacement distance of the electron from the central 
position. Providing the very basic models of the concept of refractive index and optical dispersion, 
this model has been applied by numerous authors, such as Feynman [R. P. Feynman, Lectures 
on Physics (Addison- Wesley, Massachusetts, 1963)], and Born and Wolf [M. Born and E. Wolf, 
Principles of Optics (Cambridge University Press, Cambridge, 1980)]. 

Moving on to the next level of approximation, the bound proton-electron pair may be considered 
as constituting a two-body central force problem of classical mechanics, in which one may assume 
a fixed center of mass of the system, around which the proton as well as the electron are free to 
oscillate. In this level of approximation, by introducing the concept of reduced mass for the two 
moving particles, the equations of motion for the two particles can be reduced to one equation of 
motion, for the evolution of the electric dipole moment of the system. 

The third level of approximation which may be identified is when the center of mass is allowed 
to oscillate as well, in which case an equation of motion for the center of mass appears in addition 
to the one for the evolution of the electric dipole moment. 

In each of the models, nonlinear it ies of the restoring central force field may be introduced as to 
include nonlinear interactions as well. It should be emphasized that the spring model, as now will 
be introduced, gives an identical form of the set of nonzero elements of the susceptibility tensors, 
as compared with those obtained using a quantum mechanical analysis. 

Throughout this analysis, the wavelength of the electromagnetic field will be assumed to be 
sufficiently large in order to neglect any spatial variations of the fields over the spatial extent of 
the oscillator system. In this model, the central force field is modelled by a mechanical spring force 
with spring constant fee, as shown schematically in Fig. 5, and the all-classical Newton's equations 
of motion for the electron and nucleus are 

^e-g^ = -eE{t) - ko{Xe - Xn) + h {x^ - X^)"^^ , 
optical spring 

mn-^^ = +eE{t) + ko{Xe - Xn) - ki{Xe - X^)\ 
optical spring 

corresponding to a system of two particles connected by a spring with spring "constant" 

Q^(spring) ^^(spring) 

k = -^r^ r = ^7-^^ r = ko - 2ki{xe - x„). 

d{Xe - Xn) d[Xe - Xn) 

By introducing the reduced mass^^ = ni^m^/ {rrie + m„) of the system, the equation of motion 
for the electric dipole moment p = — e(xe — a:„) is then obtained as 

d'^p ko ki o , 
ot-' rur enir nir 

This inhomogeneous nonlinear ordinary differential equation for the electric dipole moment is the 

primary interest in the discussion that now is to follow. 

The electric dipole moment of the anharmonic oscillator is now expressed in terms of a pertur- 
bation series as 

p{t) = (i) + it) + (t) + p(3) {t) + ..., 
(xE(t) (xE^{t) (xE^{t) 

where each term in the series is proportional to the applied electrical field strength to the power 
as indicated in the superscript of repective term, and formulate the system of n -f- 1 equations 
for p^''\ k = 0, 1,2, ... ,n, that define the time evolution of the electric dipole. By inserting the 



^2 Herbert Goldstein, Classical Mechanics, 2nd ed. (Addison- Wesley, Massachusetts, 1980). 
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perturbation series into Eq. (6), we hence have the equation 

q2 {0) Q2 (1) Q2 (2) ^2 (3) 
1 =- 1 1 1- 

df^ dt^ dt^ 

+ !^piO) h^pH) + hLp{2) b.p(3) + _ _ ^ 

rrir rrir rrij m-c 

+ + + p(^) + . . + P^'^ + P^'^ + •••) = —m- 

erriT 

Since this equation is to hold for an arbitrary electric field E{t) , that is to say, at least within the 
limits of the validity of the perturbation analysis, each set of terms with equal power dependence 
of the electric field must individually satisfy the relation. By sorting out the various powers and 
identifying terms in the left and right hand sides of the equation, we arrive at the system of 
equations 



^2^(0) 


+ 


^p(o) 
nil 


+ 


iLp(o)2 = 0, 
errir 






+ 


^0 (1) 

nir 


+ 


errir 


p2 


02p(2) 


+ 


_^p(2) 

rrir 


+ 


^(2pWp(2) 
errir 




9*2 


+ 




+ 


^(2p(0)p(3) 
errir 





where wc kept terms with powers of the electric field up to and including order three. At a first 
glance, this system seem to suggest that only the first order of the perturbation series depends on 
the applied electric field of the light; however, taking a closer look at the system, one can easily 
verify that all orders of the dipolc moment is coupled directly to the lower order terms. The system 
of equations for p^'^^ can now be solved for k = 0, 1,2, . . ., in that order, to successively provide 
the basis of solutions for higher and higher order terms, until reaching some k = n after which we 
may safely neglect the reamaining terms, hence providing an approximate solution. 

The zeroth order term in the perturbation series is decribed by a nonlinear ordinary differential 
equation of order two, a so-calles Riccati equation, which analytically can be solved exactly, either 
by directly applying the theory of Jacobian elliptic integrals of by applying the Riccati transor- 
mation.^'* However, by considering a system starting from rest, at a state of equilibrium, we can 
immediately draw the conclusion that p^'^^(i) must be identically zero for all times *. This, of 
course, only holds for this particular model; in many molecular systems, such in water, a perma- 
nent static dipole moment is present, something that is left out in this particular spring model 
of ours. (Not to be confused with the static polarization induced by the electric field, which by 
definition of the terms in the perturbation series is included in higher order terms, depending on 
the power of the electric field.) 

The first order term in the perturbation series is the first and only one with an explicit de- 
pendence of the electric field of the light. Since the zeroth order perturbation term is zero, the 
differential equation for the first order term is linear, which simplifies the calculus. However, since 
it is an inhomogeneous differential equation, we must generally look for a total solution to the equa- 
tion as a sum of a homogeneous solution (with zero right hand side) and a particular solution (with 
the electric field in the right hand side present). The homogeneous solution, which will contain 
two constants of integration (since we are considering second-order ordinary differential equation) 
will though only give the part of the solution which depend on initial conditions, that is to say. 



It should though be emphasized that in the limit n ^ oo, the described theory still is an exact 
description of the motion of the electric dipole moment within this model of interaction between 
light and matter. 

For examples of the application of the Riccati transformation, see Zwillinger, Handbook of 
Differential Equations, 2nd ed. (Academic Press, Boston, 1992). 
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in this case a harmonic natural oscillation of the spring system which in the presence of damping 
terms rapidly would decrease to zero. This implies that in order to find steady-state solutions, in 
which the oscillation of the dipole moment directly follows the oscillation of the electric field of the 
light, wc may directly start looking for the particular solution. For a time harmonic electric field, 
here taken as 

E{t) = E^sm{ujt), 

the particular solution for the first order term is after some straightforward algebra given as^^ 

For a material consisting of N dipoles per unit volume, and by following the conventions for the 
linear electric susceptibility in SI units, this corresponds to a first order electric polarization density 
of the form 

P«(t) = P^i)sin(a.'f) 

with the first order (linear) electric susceptibility given as 

(1)/ X (1)/ \ ^ (e'^/m,) 



where the resonance frequency ft^ = ko /nir was introduced. The Lorenzian shape of the frequency 
dependence is shown in Fig. 6. 



For the sake of self consistency, the general solution for the first order term is given as 

= Acos{{ko/mry/H) + Bsm{{kolm,f'h) + — ^^^^^ sin(a;i), 

where A and B are constants of integration, determined by initial conditions. 
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First order susceptibility 




Figure 6. Lorenzian shape of the Hnear susceptibility x^^\~^\^)- 

Continuing with the second order perturbation term, some straightforward algebra gives that 
the particular solution for the second order term of the electric dipole moment becomes 



1 



fcie^ 



fcie^ 



1 



In terms of the polarization density of the medium, still with A'' dipoles per unit volume and 
following the conventions in regular SI units, this can be written as 



p(2)(t) = p(0)+p(0)gi^(2a;f) 

= £oX^^^ (0; -uj)E^E^ + eQX^^\-2w; oj, u)E^E^ sm{2w) 

V V ' ^ V ' 

DC polarization second harmonic polarization 

with the second order (quadratic) electric susceptibility given as 



Nkie^ 
Nkie^ 



1 



1 



£0 m? (n2_a;2)(02_4^2)- 



Prom this we may notice that for one-photon resonances, the nonlinearities are enhanced whenever 
w Ri O or 2w « f2, for the induced DC as well as the second harmonic polarization density. 

The explicit frequency dependencies of the susceptibilities x^^H~2w;a;,a;) and x^^H^i ^'^) 
are shown in Figs. 7 and 8. 
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Second order susceptibility (SHG) 




0.5 1 1.5 

Figure 7. Lorenzian shape of the hnear susceptibiUty x^'^\—2uj;uj,uj) (SHG). 



Second order susceptibiiity (DC rectification) 




0.5 1 1.5 



Figure 8. Lorenzian shape of the linear susceptibility x^^HOi'^' ~'^) (-DC). 

A well known fact in electromagnetic theory is that an electric dipole that oscillates at a certain 
angular frequency, say at 2a;, also emits electromagnetic radiation at this frequency. In particular, 
this implies that the term described by the susceptibility x^^^(~2w; will generate light at 
twice the angular frequency of the light, hence generating a second harmonic light wave. 
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